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We develop a theory of normal-metal - superconductor (NS) and superconductor - superconduc- 
tor (SS) tunnelling in "bosonic" superconductors with strong attractive correlations taking into 
account coherence effects in single-particle excitation spectrum and disorder. The theory accounts 
for the existence of two energy scales, their temperature and doping dependencies, asymmetry and 
inhomogeneity of tunnelling spectra of underdoped cuprate superconductors. 

PACS numbers: 71.38.-k, 74.40.+k, 72.15.Jf, 74.72.-h, 74.25.Fy 



S3 ' 
O 1 
CD , 
i 

=5 



T3 

C 

O 

> 
in 



o 



X 



Soon after the discovery of high-T c superconductiv- 
ity Qi a number of tunnelling, photoemission, optical, 
nuclear spin relaxation and electron-energy-loss spectro- 
scopies discovered an anomalous large gap in cuprate 
superconductors existing well above the superconduct- 
ing critical temperature, T c . The gap, now known as 
the pseudogap, was originally assigned [2] to the binding 
energy of real-space preformed hole pairs - small bipo- 
larons - bound by a strong electron-phonon interaction 
(EPI). Since then, alternative explanations of the pseu- 
dogap have been proposed, including preformed Cooper 
pairs [H, inhomogeneous charge distributions containing 
hole-rich and hole-poor domains 0], or competing quan- 
tum phase transitions [5j. 

Present-day scanning tunnelling (STS) 0, 0, S|, in- 
trinsic tunnelling [9( and angle-resolved photoemission 
(ARPES) [5J] spectroscopies have offered a tremendous 
advance into the understanding of the pseudogap phe- 
nomenon in cuprates and some related compounds. Both 
extrinsic (see @, Q and references therein) and intrin- 
sic @ tunnelling as well as high-resolution ARPES [|[ 
have found another energy scale, reminiscent of a BCS- 
like "superconducting" gap that opens at T c accompa- 
nied by the appearance of Bogoliubov-like quasi-particles 
Q around the node. Earlier experiments with a time- 
resolved pump-probe demonstrated two distinct gaps, 
one a temperature independent pseudogap and the other 
a BCS-like gap Also, Andreev reflection experi- 

ments revealed a much smaller gap edge than the bias at 
the tunnelling conductance maxima in a few underdoped 
cuprates (llj . Another remarkable observation is the spa- 
tial nanoscale inhomogeneity of the pseudogap observed 
with STS [(| 0, HI and presumably related to an unavoid- 
able disorder in doped cuprates, Fig. la. Essentially, the 
doping and magnetic field dependence of the supercon- 
ducting gap compared with the pseudogap and their dif- 
ferent real space profiles have prompted an opinion that 
the pseudogap is detrimental to superconductivity and 
connected to a quantum critical point rather than to pre- 
formed Cooper pairs 0] . Nevertheless without a detailed 
microscopic theory that can describe highly unusual tun- 
nelling and ARPES spectra, the relationship between the 
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FIG. 1: (Colour online) Ratio, R = 7„ s (-100)//„ s (100), 
of the negative-bias NS tunnelling conductance to the pos- 
itive bias conductance fl2l ] integrated from meV to ^100 
meV, respectively, for a few cuprate superconductors in a 
wide range of atomic hole density, p. Inset (a): atomi- 
cally resolved asymmetric STS spectra of Lai.88Sro.i2Cu04 
at 4.2K acquired at different points of the scan area Q|; (b): 
momentum-integrated photoemission showing no vHs of DOS 

EH. 



pseudogap and the superconducting gap remains a mys- 
tery Q . It has become increasingly likely that NS and SS 
tunnelling spectra of underdoped cuprates do not agree 
with the simplest BCS spectra neither s-wave nor d-wave. 
Apart from the almost temperature independent pseudo- 
gap, A p with 2A p / 'ksTc often many times larger than 
the BCS ratio (ss 3.5), there is an asymmetry in NS 
tunnelling, Fig.l. The integrated conductance for the 
negative bias is larger than for the positive bias in many 
samples. The van Hove singularity (vHs) of the density 
of states (DOS) is ruled out as a possible origin of the 
asymmetry since it is absent in the momentum-integrated 
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photoemission, Fig. lb. The asymmetry is expected for 
conventional semiconductors or Mott-Hubbard insulators 
(Ins(— oo)/I na (+00) = (1 — p)/2p), but neither of them 
account for its magnitude, Fig.l, if disorder and matrix 
elements are not considered. Here we develop the theory 
of NS and SS tunnelling in bosonic superconductors with 
strong attractive correlations by taking into account 
disorder and coherence effects in a single- particle excita- 
tion spectrum. Our theory accounts for peculiarities in 
extrinsic and intrinsic tunnelling in underdoped cuprate 
superconductors. 

Recent Monte Carlo calculations show that high-T c su- 
perconductivity cannot be explained by the simplest re- 
pulsive Hubbard model [H[ , so that one has to extend the 
model to get some superconducting order [T3|. On the 
other hand even a moderate EPI significantly increases 
the superconducting condensation energy [17| stabiliz- 
ing mobile small bipolarons [HI, 0], as anticipated for 
strongly correlated electrons in highly polarizable ionic 



lattices [1J]. Real-space pairs, whatever the pairing in- 



teraction is, can be described as a charged Bose liquid on 
a lattice, if the carrier density is relatively small [14( . The 
superfluid state of such a liquid is the true Bose-Einstein 
condensate (BEC), rather than a coherent state of over- 
lapping Cooper pairs. Single-particle excitations of the 
liquid arc thermally excited single polarons propagating 
in a doped insulator band or are localised by impurities. 
Different from the BCS case, their negative chemical po- 
tential, fi, is found outside the band by about half of 
the bipolaron binding energy, A p , both in the supercon- 
ducting and normal states Here, in the supercon- 
ducting state (T<T C ), following Ref.[20|] we take into ac- 
count that polarons interact with the condensate via the 
same potential that binds the carriers, so that the single- 
particle Hamiltonian in the Bogoliubov approximation is 



p = E 



LrP'rPr + ( -A,,,/^//,. • H.C 



(1) 



where — E u — fi, E v is the normal-state single-polaron 
energy spectrum in the crystal field and disorder po- 
tentials renormalised by EPI and spin-fluctuations, and 
A c „ = — A cP is the coherent potential proportional to 
the square root of the condensate density, A c oc ^Jn c (T). 

The operators p\, and p\, create a polaron in the single- 
particle quantum state v and in the time-reversed state 9, 
respectively. As in the BCS case the single quasi-particle 
energy spectrum, e„, is found using the Bogoliubov trans- 
formation, p v = u u a u + v w (3l, p 9 = u v (3 w - v u al, 

e v = [e v + Ay 1/2 , with ul vl = (l±&/e„)/2. This 
spectrum is different from the BCS quasi-particles be- 
cause the chemical potential is negative with respect to 
the bottom of the single-particle band, /i = — A p . A 
single-particle gap, A, is defined as the minimum of 
e„. Without disorder, for a point-like pairing poten- 
tial with the s-wave coherent gap, A c k ~ A c , one has 
H3 A(T) = [Ap + A C (T) 2 ] 1/2 . The full gap varies with 

1/2 

temperature from A(0) = [A 2 + A c (0) 2 ] at zero tem- 
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FIG. 2: (Colour online) LDA+GTB energy band structure of 
underdoped cuprates with the impurity localised states [22I ] 
shown as horizontal lines in (a); NS model densities of states 
(b) showing the bandtail in the bosonic superconductor. 



perature down to the temperature independent A = A p 
above T c , which qualitatively describes some earlier and 
more recent [§] observations including the Andreev re- 
flection in cuprates (see [2(| and references therein) . The 
NS and SS tunnelling transitions are described with the 
tunnelling Hamiltonians, H ns and H ss respectively, [2l| . 
which are perturbations: 

H ns = pJ2pI c - + BN- 1 / 2 b^ppcu + h.c, 

vv' vv'i]' 

h ss = pY,pIpu + bn- 1 / 2 x 



E + E b l'P»pv + h - c - 



(2) 



n' 



Here c„, p u > and b 7 y annihilate a carrier in the normal 
metal, a single polaron and a composed boson in the 
bosonic superconductor respectively, N is the number of 
unit cells. Generally the tunnelling matrix elements with 
(P) and without (P) involvement of the composed boson 
are different, B > P, because the presence of an addi- 
tional hole lowers the tunnelling barrier for an injection 
of the electron [2l| . Applying the Bogoliubov transfor- 
mation to Eqs.Q and the standard perturbation theory 
yields the following current-voltage characteristics: 

Ins{V) = ^-^f- Y> 2 AFv ~ fvm* + eV-e u ,) + 



v 2 v ,(F u + f u ,-l)5(Z v + eV + e v ,)} + 



2neB 2 



E 



[ul,[F v f v , - {x/2){\ -F v - +eV + e v >)+ 

v 2 v ,[F v {l - /„,) + (x/2)(F u - /„0W& + eV - e„)] , 

(3) 

where F v = l/[exp(^/fc S T) + 1], /„, 
l/[exp(e l/ '/fesT) + 1] are distribution functions of 
carriers in the normal metal and single quasi-particles 



3 



(4) 



respectively, and x/2 is the atomic density of composed 
bosons in the superconductor, and 

5(e u + eV -e v ,) + ulvl,{f v + f v , -l)x 
[S(e v + eV + e„,)- 5(e v - eV + €„>)]]+ 

2 -^f- Yy v ul,{{l -f v - f v ,)x/2 - /„/„,)+ 

vlv 2 A{l- f v - f u ,)x/2 + (1 - /„)(1 - f v ,))] 
x [8{e v - eV + e v >) - 5{e v + eV + e v ,)} + 

2u 2 v v 2 v ,[{f vl - f v )x/2- f u (l- f v ,)] 
x [8{e v + eV - e v >) - 6(e„ - eV - e v >)], 

where V is the voltage drop across the junction. For 
more transparency we neglect the boson energy disper- 
sion in Eqs.© 0J, assuming that bosons are sufficiently 
heavy, so their bandwidth is relatively small. Here 
we adopt the "LDA+GTB" band structure with im- 
purity bandtails near (tt/2,it/2) of the Brillouin zone, 
Fig. 2a, which explains the charge-transfer gap, E ct , the 
nodes and sharp "quasi-particle" peaks, and the high- 
energy "waterfall" seen in ARPES [221 ]. The chemical 
potential is found in the single-particle bandtail within 
the charge-transfer gap at the bipolaron mobility edge, 
Fig. 2b, in agreement with the SNS tunnelling experi- 
ments [23] - Such a band structure explains an insulating- 
like low temperature normal-state resistivity as well as 
many other unusual properties of underdoped cuprates 
[14j . If the characteristic bandtail width of DOS, T, 
is sufficiently large compared with the coherent gap, 
r > A c „, one can factorize the quasi-particle DOS as 
P(E) = £„ S(E - e v ) « [p n (E) + p n {-E)] Pa {E) for any 
symmetry of the coherent gap. Here p n (E) is the nor- 
mal state DOS of the dop ed insulat or with the band-tail, 
Fig. 2a, and p s {E) = Ej \J E 2 — A 2 for the s-wave gap, or 
p a (E) = (2/TT)[e(l-E/A )EK(E/A )/A + e(E/A Q - 
1)K(A /E)] for a d-wave gap, A cv — A cos(2^) [24[ 
(K(x) is the complete elliptic integral and 4> is an angle 
along the constant energy contour). One can neglect the 
energy dependence of the normal metal DOS. Then dif- 
ferentiating Eq. ([3|) over the voltage and integrating yields 
the NS tunnelling conductance a ns = dI ns /dV at zero 
temperature for the d-wave case, 
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FIG. 3: (Colour online) Theoretical NS conductance, Eq.((5]l, 
for A = T, A p = 2.7F and B = 2.65P. The superconducting 
gap and the pseudogap are shown with triangles and arrows, 
respectively. Inset shows a representative STS spectrum of 
I^-zSr^CuCU with x = 0.12 at 4.2K 0. 



and the model parameters can be made without affect- 
ing our conclusions as long as the characteristic features 
are reflected in this choice. Eq.© captures all unusual 
signatures of the experimental tunnelling conductance in 
underdoped cuprates, such as the low energy coherent 
gap, the high-energy pseudogap, and the asymmetry. In 
the case of atomically resolved STS one should replace 
the averaged DOS p n (E) in Eq.© with a local bandtail 
DOS p n (E,r), which depends on different points of the 
scan area r due to a nonuniform dopant distribution. As 
a result the pseudogap shows nanoscale inhomogeneity, 
while the low-energy coherent gap is spatially uniform, 
as observed [1], Fig. la. Increasing doping level tends to 
diminish the bipolaron binding energy, A p , since the pair- 
ing potential becomes weaker due to a partial screening 
of EPI with low-frequency phonons [25[. However, the 
coherent gap, A c , which is the product of the pairing po- 
tential and the square root of the carrier density [2C| , can 
remain about a constant or even increase with doping, as 
also observed Q. 

In the case of the SS tunnelling we use Eq.© to ad- 
dress two unusual observations: a gapped conductance 
near and above T c and a negative excess resistance be- 
low T c [i| [26| . Eq. © is grossly simplified in the normal 
state, where A cv — 0, 

iss(v) - YJJ" - + eV ~ 



o ns oc A + p s (\eV\)[p n (-eV) + p n (eV)} + 
A-[l - 2cos' 1 (|ey|/A )O(l - \eV\/A )/n] 
x [ Pn (-eV) - Pn (eV)}, (5) 

where A ± = 1 ± B 2 [9(— eV) + x/2] /P 2 , and B(E) is 
the Heaviside step function. The theoretical conduc- 
tance, Eq.©, calculated with a model normal state DOS, 
p n {E)/pb = {l+tanh[(_B-Ap)/r]}/2, and the d-wave su- 
perconducting DOS, p s (E), is shown in Fig. 3. Our model 
p n {E) reflects the characteristic energy dependence of 
DOS in disordered doped insulators, which is a constant 
Pb above the two-dimensional band edge and an expo- 
nent deep in the tail. Any particular choice of p n {E) 
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x [8(£ v -eV + &)- +eV + 



(6) 



Near and above the transition but sufficiently below the 
pseudogap temperature T* = A p /kB > T > T c , and if 
the voltage is high enough, eV > fcgT, one can neglect 
temperature effects in Eq.© and approximate f v with 
the step function, j v = Q(— £, v ). So using the model 
normal state DOS yields 
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FIG. 4: (Colour online) Approximate normal state tunnelling 
conductance of bosonic superconductor (solid line) with V = 
3.2 meV and A p = 16 meV compared with the experimen- 
tal conductance [9J (symbols) in mesas of Bi2Si'2CaCu208+a 
(T C =95K) at T = 87K. 



B 2 (l + x/2) 1 + a 2 b 2 | B 2 xa 2 a 2 + b 2 

P 2 {a 2 b A - 1) n a(l + b 2 ) + 2P 2 (a 2 - 6 4 ) n a(l + 6 2 ) ' 

(7) 

where a = exp(|eV~| /T) and & = exp(A p /r). When 
b ^> 1 and x is not too small, the first two terms 
on the right hand-side of this equation are negligible. 
Hence the tunnelling matrix elements and doping have 
little effect on the shape of the current-voltage depen- 
dence. At sufficiently high voltages eV > fc^T the 
conductance (from a(V) — dI ss /dV with Eq.([7j)) ac- 
counts for the gaped conductance in underdoped mesas 
of Bi2Sr2CaCu2 0g+5 near and above T c , as shown in 
Fig 4. The finite temperature neglected in Eq.([7]) ac- 
counts for some excess experimental conductance at low 
voltages in Fig. 4 compared with the theoretical conduc- 
tance. The negative excess resistance below T c [26[ 



can be explained by expanding Eqs. (|4I6[ ) in powers of 
eV giving a zero bias conductance. For low temper- 
atures in the superconducting state this is er s (0) oc 
Io° dep s (e) 2 cosh(e/2fc B T)~ 2 , and in the normal 
state (J„(0) oc T' 1 d^p n {^) 2 cosh^/2k B T)- 2 . Es- 
timating these integrals yields, respectively a s (0) oc 
T -1 exp(— Ac/ksT) for the s-wave coherent gap, or 
cr s (0) oc T 2 for the d-wave gap, and er n (0) oc 
T _1 exp(— T*/T). The latter expression is in excellent 
agreement with the temperature dependence of the mesa 
tunnelling conductance above T c [26| (see also Ref. [25j). 
Extrapolating this expression to temperatures below T c 
yields the resistance ratio R s /R n oc exp[(A c /fc B — T*) /T] 
(s-wave) or R s /R n oc exp(-T* /T)/T 2 (d-wave). Hence 
in underdoped cuprates, where T* > A c /fcs, the zero- 
bias tunnelling resistance at temperatures below T c is 
smaller than the normal state resistance extrapolated 
from above T c to the same temperatures (i.e. the nega- 
tive excess resistance), as observed [26j |. 

In summary, we have developed the theory of tun- 
nelling in bosonic superconductors by taking into ac- 
count coherence effects in the single-quasi-particle energy 
spectrum, disorder and the realistic band structure of 
doped insulators. The theory accounts for the existence 
of two energy scales in the current-voltage NS and SS 
tunnelling characteristics, their temperature and doping 
dependence, and for the asymmetry and inhomogeneity 
of NS tunnelling spectra of underdoped cuprate super- 
conductors. 
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